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ABSTRACT 

We study statistical properties of long gamma-ray bursts (GRBs) produced by the 
collapsing cores of WR stars in binary systems. Fast rotation of the cores enables a 
two-stage collapse scenario, implying the formation of a spinar-like object. A burst 
produced by such a collapse consists of two pulses, whose energy budget is enough to 
explain observed GRBs. We calculate models of spinar evolution using results from 
a population synthesis of binary systems (done by the 'Scenario Machine') as initial 
parameters for the rotating massive cores. Among the resulting bursts, events with the 
weaker first peak, namely, precursor, are identified, and the precursor-to-main-pulse 
time separations fully agree with the range of the observed values. The calculated 
fraction of long GRBs with precursor (about 10 per cent of the total number of long 
GRBs) and the durations of the main pulses are also consistent with observations. 
Precursors with lead times greater by up to one order of magnitude than those observed 
so far are expected to be about twice less numerous. Independently of a GRB model 
assumed, we predict the existence of precursors that arrive up to > 10'^ s in advance 
of the main events of GRBs. 
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1 INTRODUCTION 

Gravitational collapse is believed to be the underlying mech- 
anism for the most energetic events observed in the Universe: 
GRBs and supernovae. While it is commonly accepted that 
the remnant of such events is a black hole or a neutron 
star, the details of the process are uncertain. In relation to 
GRBs, we investigate here the collapse of a fast rotating 
magnetized object, which can be understood in terms of the 
'spinar paradigm'. We define spinar as a critically-fast rotat- 
ing magnetized relativistic object, whose quasi-equilibrium 
is maintained by the balance of centrifugal and gravitational 
forces. The evolution of a spinar is determined by its mag- 
netic field. A benefit of the spinar model is that it describes 
transparently and in a simple way the main features of a 
real collapse. 

The properties of rotating magnetized objects were first 
investigated to understand the mechanisms of active galac- 
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tic nuclei by, e.g., Hoyle & Fowler (1963); Ozernoi (1966); 
Morrison (1969); Woltjer (1971); Bisnovatyi-Kogan & Blin- 
nikov (1972); Ozernoy & Usov (1973), of pulsars by Gunn 
& Ostriker (1969), and of supernova explosions by LeBlanc 
& Wilson (1970); Bisnovatyi-Kogan (1971). A rotation- 
supported 'cold' configuration with magnetic field received 
the name 'spinar' (see early reviews by Morrison & Cava- 
liere 1971; Ginzburg & Ozernoi 1977). Stellar mass spinars 
were suggested by Lipunov (1983, 1987). In the works by 
Lipunova (1997) and Lipunova & Lipunov (1998) a burst of 
electromagnetic radiation produced during the collapse of a 
spinar was studied, and a spinar mechanism for GRBs was 
first suggested. 

As Lipunov & Gorbovskoy (2007) point out, there 
should be energy release in a process of spinar formation 
as well. This approach enables one to consider a two-stage 
scenario of a collapse. At the first stage, a spinar forms from 
a collapsing rotating body when centrifugal forces halt con- 
traction. The effective dimensionless Kerr parameter of the 
spinar is greater than unity. At the second stage, the angu- 
lar momentum is carried away, and the spinar evolves to a 



2 G. V. Lipunova et al. 



limiting Kerr black hole or a neutron star, depending on its 
mass. 

Lipunov & Gorbovskoy (2008) develop a ID model of 
the magneto-rotational collapse of a spinar, which includes 
all principle rclativistic effects on the dynamics and the mag- 
netic field, along with the pressure of nuclear matter and 
neutrino cooling. A variety of burst patterns is obtained, 
generally a combination of two peaks. It is shown that the 
spinar paradigm agrees with the basic observed GRB prop- 
erties. 

Potential progenitors of spinars are rotating WR stars 
without H and He envelops, which are already considered as 
possible progenitors of GRBs (for a review, sec Woosloy & 
Bloom 2006). The spinar mechanism requires the presence of 
a high angular momentum in the WR core at the start of the 
collapse. A direct collapse to a black hole is impossible if the 
rotating core has an eilective dimensionless Kerr parameter 
greater than unity: 



where I is the moment of inertia of the core, n is the angular 

velocity, Mc is the core mass. Equation (1) corresponds to 
the following condition on the specific angular momentum: 

in , , ,„15 Mc 2-1 

— > 4.4 X 10^^ — ^ cm s \ 
Mc Mq 

Fast rotation can be a result of the evolution of a rotating 
single massive star or a star in a close binary system (see, 
e.g., van den Heuvel & Yoon 2007). 

Different scenarios for the collapse of a WR core arc pos- 
sible depending on the unknown properties of the collapsing 
core, among them the quantity and distribution of the an- 
gular momentum within the core. These scenarios have so 
far included: a highly magnetized rotating neutron star, a 
black hole surrounded by an accretion disc (the 'coUapsar' 
model), a hypermassive rotating neutron star with an accre- 
tion disc (see references in Woosley & Bloom 2006). 

In this regard, the spinar mechanism is a natural com- 
plement to the range of potential GRB producers. It is im- 
portant to note that in distinction from the coUapsar model 
of Woosley (1993), in the spinar paradigm a GRB begins 
with a spinar formation and continues with its collapse to 
a hmiting black hole. In the coUapsar model, we first have 
the formation of a black hole, and after that a GRB devel- 
ops, powered from an accreting disc-like envelope and by the 
Blandford-Znajek mechanism. It is essential for the spinar 
that the central part of a core with mass ~ 2— 3Mq has large 
angular momentum (effective Kerr parameter > 1). In con- 
trast, the coUapsar model requires that there is less angular 
momentum in the centre and an excess at the periphery. 

There are calculations of the structure of WR cores sup- 
porting a hypothesis that may be there is too much angular 
momentum in the centre. For example, the results of Hirschi 
et al. (2005); Yoon et al. (2006) indicate that the inner part 
of a WR core with mass ~ 2.5 M© is characterized by the 
effective dimensionless Kerr parameter not less (or not sig- 
nificantly less) than an effective Kerr parameter of the whole 
core. Consequently, if a core has the effective Kerr param- 
eter considerably greater than unity and cannot undergo a 
direct collapse to a black hole, then the inner part is equally 
prohibited to do so. The results of such models, usually in- 



tended for single stars, are generally uncertain as they are 
strongly dependent on the physics assumed. 

In the present work we consider only massive stars col- 
lapsing in binary systems. In a binary system, tidal interac- 
tion and synchronization load to the fast rotation of a mas- 
sive prc-coUapse core (Tutukov & Yungelson 1973). Binary 
systems as GRB sites were considered by, e.g., Woosley 
(1993); Paczynski (1998); Brown et al. (2000); Postnov & 
Cherepashchuk (2001); Tutukov & Cherepashchuk (2003, 
2004); Izzard et al. (2004); Podsiadlowski et al. (2004); 
Petrovic et al. (2005); Bogomazov et al. (2007); van den 
Heuvel & Yoon (2007). 

In previous work (Lipunov & Gorbovskoy 2008), the 
dependence of qualitative features of GRBs on the initial 
parameters of collapsing cores (their angular momentum, 
magnetic flux, and mass) have been studied by calculating 
numerous models on a uniform grid of parameter values. 
It is clear, however, that the initial parameters are not dis- 
tributed uniformly, and their values are defined by the course 
of previous evolution. 

To derive the initial distribution on the efi^ective Ker 
parameter, we use a population synthesis of binary stars 
carried out by the Scenario Machine (Lipunov et al. 1996a, 
2007). We perform population synthesis for different binary 
evolution parameters and analyse the resulting set of GRBs, 
which can be classified into different kinds of events. 

Below we focus particularly on the precursor phe- 
nomenon (Koshut et al. 1995; Lazzati 2005; Burlon et al. 
2008; Wang & Meszaros 2007, and references therein), re- 
garding it as a primary energy release accompanying a spinar 
formation. The 'lead time' of a precursor, namely, the sep- 
aration time from the main pulse, is an important and ro- 
bust parameter that we calculate. As we show, there should 
be GRBs with precursor occuring hunderds and thousands 
of seconds in advance. We also investigate the dependences 
between energy and temporal characteristics of two GRB 
pulses. 

In the next section we give a brief account of the rel- 
evant aspects of the population synthesis of binary stars. 
In Sect. 3 we describe calculation of spinar evolution. Re- 
sults on GRB statistics are presented in Sect. 4. Discussion 
follows in Sect. 5. 



2 POPULATION SYNTHESIS OF WR STARS 
IN BINARY SYSTEMS 

Population synthesis of binary systems is performed by the 
'Scenario Machine', described in detail by Lipunov et al. 
(1996a, 2007). Population synthesis of WR stars in binary 
systems as progenitors of long GRBs is done by Bogomazov 
et al. (2007, 2008) to investigate the rates and correlation 
with the host-galaxy morphology. Here we describe only a 
few elements of the computer code, which involve the pa- 
rameters varied in the present study. 

The minimum initial mass of a star that may even- 
tually produce a black hole is set to 25M0 (Tutukov & 
Cherepashchuk 2003).^ Its companion can have any mass be- 

^ If in the course of the evolution a star accretes matter from 
the companion, we check the maximum mass it attains after the 
accretion. 



Population synthesis of gamma-ray bursts and the spinar paradigm 3 



tween 0.1 Mq and the mass of the main star, so the mass ra- 
tio q = M2/M1 < 1. The initial masses of the main stars are 
distributed by Salpeter's law above 25Mq. As representa- 
tive cases, two initial distributions of mass ratio f{q) = g"' 
are considered; with aq — and Qg — 2. 
Two models of a stellar wind arc used: 

Wl: 'Weak stellar wind'. This scenario corresponds to 
that named 'A' in Bogomazov et al. (2007); Lipunov et al. 
(2007). On the main sequence and at the supergiant stage, 
the total mass loss does not exceed 20 per cent of the initial 
mass. During the WR-star stage, the mass loss is 30 per cent 
of the whole star mass. Weak stellar winds are expected for 
stars with low metallicity. 

W2: 'Strong stellar wind'. This scenario corresponds to 
that named 'C in Bogomazov et al. (2007); Lipunov et al. 
(2007), except that a star loses 70 per cent of its envelope at 
each evolutionary stage, not the whole envelope. Note that, 
in the original model 'C, GRBs cannot be produced because 
the strong wind results in a too large semi-major axis of the 
orbit before the collapse of a binary component (Bogomazov 
et al. 2007). 

In model Wl, the mass loss rate M for a main sequence 
star is described by the formula M = Ow Lc/Vo^ , where L 
is the star luminosity, Voo is the wind velocity at infinity, c 
is the speed of light, a„ is a free parameter. On the main 
sequence and at the supergiant stage, the mass variation 
in model Wl during any stage does not exceed a value of 
0.1 (M — Mcore), where M is the mass of the star at the 
beginning of the stage and Mcore is the mass of the stellar 
core. 

During the common- envelope stage, underwent by a 
quarter to a half of the systems, stars give their angu- 
lar momentum very effectively to the surrounding mat- 
ter and spiral towards each other. The effectiveness of 
the common-envelope stage is described by the parameter 
acE = A_Bb/Ai?orb, where AEh = iSgrav — -Bthermai is the 
binding energy of the matter lost to the envelope and AEoib 
is the decrease in the energy of gravitational interaction 
of the approaching stars. The smaller parameter acE, the 
closer components become after the CE stage. 

Further details on the wind scenarios and restrictions 
on the parameters of the binary evolution can be found 
in Lipunov et al. (1996b, 1997, 2005). In the present study, 
we take typical values of the parameters, also used in Bogo- 
mazov et al. (2007). 

To estimate the effective Kerr parameter at the start 
of collapse, we need a relation between the radius and the 
mass of a rotationally-synchronized core at the evolution- 
ary stage, after which, and until the collapse, the angu- 
lar momentum is conserved. Following previous studies, we 
believe that, at the He-burning stage, the star rotation is 
fully synchronized with its orbital motion (e.g., Tutukov & 
Cherepashchuk 2003, 2004; Izzard et al. 2004; Podsiadlowski 
et al. 2004; van den Heuvel & Yoon 2007; Belczynski et al. 
2007; Zahn 2008). We take into consideration that the car- 
bon burning is likely to be completed before the components 
of even a very close binary are re-synchronized (Masevich & 
Tutukov 1988; Zahn 2008). Thus, during and after the C- 
burning stage, the period of the axial rotation of the core is 
less than the orbital period. 

One can estimate the radius of the core by the end of 



Table 1. Parameters for scenarios of the population synthesis 
and the resulting ratio / of the number of the cores with oq > 1 

to the total number of the cores. 



Population 


Wind model 




"CE 


aq 


ki 


/ 


1 


Wl 


0.3 


0.5 





0.4 


0.10 


la 


Wl 


0.3 


0.5 





1.0 


0.21 


lb 


Wl 


0.3 


0.5 





0.1 


0.02 


2 


Wl 


0.1 


0.5 





0.4 


0.10 


3 


Wl 


0.7 


0.5 





0.4 


0.09 


4 


Wl 


0.3 


0.5 


2 


0.4 


0.05 


5 


Wl 


0.3 


1.0 





0.4 
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helium burning from the mass and temperature using the 
virial theorem: 

6kT ' 

where -Rc and Mc are the core radius and mass; T is the 
temperature of carbon burning, which is about 6 X 10* K 
(Masevich & Tutukov 1988); G, the gravitational constant; 
IX, the average number of nucleons for a particle (equal to 
15); rup, the proton mass; k, the Boltzmann constant. Here 
we use the fact that CO-core of a massive star is not degen- 
erate. 

The effective Kerr parameter is calculated following 
definition (1). The moment of inertia can be written as 
I — ki Mc Rc , where a dimensionless parameter ki likely 
lies within ~ 0.1 — 0.4 {ki = 0.4 for a uniform spherical 
body; kj = 0.1 for polytrope spheres (Zeldovich et al. 1981, 
chap. 2.2)). 

We perform population synthesis for different values of 
the evolutionary parameters a„, acE, aq, and different wind 
models (Table 1). Fig. 1 presents the resulting distributions 
of the CO-cores of WR-stars against the value of the effective 
Kerr parameter. 

The last column in Table 1 presents the resulting ratio 
of the number of the cores with the effective Kerr parameter 
greater than unity to the total number of the cores. Though 
the absolute rate of GRBs is not the subject of our study, 
and it cannot be inferred here, we note that the relative 
frequencies of GRB production in different binary evolution 
scenarios are reflected by the numbers in the last column 
of Table 1. An influence of the possible magnetic coupling 
between the core and the slower rotating envelope, which 
might brake core rotation up to several times after the end 
of helium burning (see, e.g., van den Heuvel & Yoon 2007), 
is illustrated by population 'lb'. For 'lb', we decrease the 
coefficient of inertia ki; effectively, this describes some loss 
of the angular momentum of the core and leads to lower 
number of fast-rotating cores. Note that our numerical esti- 
mates agree with those of Belczynski et al. (2007) who study 
the evolution of Population HI binaries^. 

^ See tabic 2 of their work, model ModOS. This model is calcu- 
lated with the zero magnetic coupling between the core and the 
envelope, for the flat initial mass ratio distribution. The fraction 
of the cores with the effective Kerr parameter of the inner 7 Mq 
greater than unity is about 0.1 — compare columns rji and 773 — 
just what we observe in our Table 1 for weak wind models '1', 
'2', '3', and '5' with fe/ of the same order as in Belczynski et al. 
(2007), who use for the exposed He-cores ki = 0.2 — 0.3. 
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Figure 1. Probability density (a) and cumulative distribution 
(b) of the effective Kerr parameter (oq > 1) of the CO-cores of 
WR stars in binary systems with different evolutionary parame- 
ters and fcj listed in Table 1. Numbers on the vortical axis are 
normalized to the total number of the cores with ao > 1, which 
is different for each population. 



For the sets of parameters, used in populations '1', '4', 
and '5', Bogomazov et al. (2007) performed population syn- 
thesis and estimated absolute rates of collapse events with 
orbital periods less than ~ 1 day. This corresponds to high 
values of the effective Kerr parameter according with 



ao 



0.3 



Mc Iday ki 



0.4 



which follows from equations (1) and (2). Close binaries con- 
taining a black hole, a main-sequence star, or a nondegener- 
ate star filling its Roche lobe contribute to this result with 
different rates. The rates obtained by them agree with the 
observed total GRB rate derived implying a solid angle cor- 
rection for GRB emission (Podsiadlowski et al. 2004). 



3 CALCULATION OF SPINAR MODELS 

3.1 A two-stage collapse of a spinar 

Here we describe the principles of the spinar paradigm (for 
more details, see Lipunov & Gorbovskoy 2007, 2008). Con- 



sider the collapse of a rotating object. When there is no suffi- 
cient internal pressure to balance the gravitational force, the 
body will contract until it encounters a centrifugal barrier 



at a radius Tisp estimated from the relation: 
2 „ GM 



(3) 



where uj is the angular velocity. At this point we say that 
a spinar forms. Actual dynamical behavior of the configu- 
ration can be very complex at this transition; however, the 
spinar model results in the correct amount of the energy re- 
leased in the process - half of the gravitational energy. From 
equation (1) one gets: 



P _ 2GM 

-Tlsp — do ^ 



(4) 



This radius is to be considered as a characteristic size in the 
equatorial plane'^. The effective Kerr parameter is virtually 
constant during the first stage. 

The evolution of a spinar proceeds as its angular mo- 
mentum decreases because of magnetic and viscous forces 
producing a braking torque. A spinar can be characterized 
by an average magnetic field which is represented by the 
magnetic dipole moment /i. The braking torque on a spinar 
can be expressed by a general formula: 



d (MR^uj) 
At 



(5) 



(see also Lipunov 1992, chapter 5), where R is the spinar 
radius and kt is a dimensionless parameter of order of unity. 
Such an approach describes a maximally effective mecha- 
nism for spin-down. Note that the spinar's rotation speeds 
up with decreasing spin. 

The magnetic dipole moment p can be expressed using 
parameter Qm = -Bq -Ro/(6 G Mq) - the ratio of the magnetic 
energy to the gravitational energy - a measure of magneti- 
zation introduced by Lipunov & Gorbovskoy (2007). The 
initial value of magnetic dipole moment /lo ~ Bo R^ /2 can 
be rewritten as y/oui G Mo Ro, incorporating constants of 
order of unity into Qm- Here zero-subscripts indicate initial 
values of the parameters, and B is the dipolar strength of the 
magnetic field, which approximates the actual field configu- 
ration. In the approximation of magnetic fiux conservation, 
the magnetic parameter Qm remains constant. 

As the spinar contracts and its rotation speeds up, the 
energy-release rate increases, but as the spinar's size ap- 
proaches the size of a limiting Kerr black hole, the magnetic 
field squeezes up against the surface, and the spinar's power 
fades away. Thus the second peak is produced. To account 
for the vanishing magnetic field, we use the following ap- 
proximate expression (Ginzburg & Ozernoi 1965; Lipunova 
1997; Lipunov & Gorbovskoy 2008): 



M = MO 



C(-Ro/rg,ao) 
C(i?/rg,a) 

n(a) 



+ 



'miii(^) 

2x2 



+ in 1 - 



(6) 



^ In the spinar equations it is assumed that kj = 1. This value 
ensures self-consistency of a calculated spinar evolution when the 
spinar's size approaches GMjc?. 
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where Vg = G M/(? and ^(a;,a) ~ — (a;min/2:)'^/3 if a; ^ a;min 
and ^(a;,a) —oo it x Smin, with Xminia.) = 1 + vT~~a 
for a < 1 and a;niin(a) ~ 1 for a ^ 1. 

Lipunov & Gorbovskoy (2008) obtain the characteris- 
tic time-scale for the second stage, or the time scale of the 
angular momentum losses (c.f. 3 and 5): 



1500 



GM 



2fct On 



(7) 



The complete set of equations (see Lipunov & Gor- 
bovskoy 2008) includes all the main relativistic effects near 
the gravitational radius: dynamics in the Kerr metric, mag- 
netic field extinction, time dilation, and gravitational red- 
shift. The 'spinar engine' posesses properties commonly con- 
sidered necessary to launch jets: rotation, magnetic field, rel- 
ativistic velocities and energies. These jets can be observed 
as GRBs. The magneto-dynamical model of the spinar pro- 
vides a main-scale pattern of time evolution for a GRB, ig- 
noring its short-scale variability. We do not make any as- 
sumptions about the specific nature of the jets, but we will 
imply for self-consistency that there is an efficiency of pro- 
cessing the burst energy to the jet, which can be about 1 
per cent (see also Sect. 5.2). 

3.2 Characteristics of spinar 'light curves' 

The input parameters to the model of collapse are the effec- 
tive dimensionless Kerr parameter of the core, the magnetic 
parameter, and its mass. For example. Fig. 2 shows the dis- 
tribution of the effective Kerr parameters and masses of the 
cores of WR stars in close binary systems for population '1' 
(throughout the paper we denote a result from population 
synthesis, i.e. a set of massive cores, as a population with 
a corresponding number from the first column of Table 1). 
Only cores with ao > 1 are counted in Fig. 2. We note that 
for populations 'l'-'5' the distribution of WR core masses 
with unrestricted ao has one peak around 7 — 8 Mq . More 
massive cores have an advantage in acquiring large ao as they 
have a greater radius while synchronized (see equation (2)); 
this results in the multi-peak form of the mass distribution 
in Fig. 2. At the same time, the mass distribution itself is 
not very important for the results of the present work (it 
only shows itself in Fig. 5a as a three-branch disposition of 
points). For each pair ao and m = M/Mq, we set a value 
of Om randomly and uniformly distributed on a logarithmic 
scale with limits 10~^..10~^. 

To calculate spinar temporal evolution, we solve a sys- 
tem of differential equations (Lipunov & Gorbovskoy 2008) 
using a fourth-order Rosenbrock method (Press et al. 2002, 
Chap. 16) with tolerance 10~^. We start integration of equa- 
tions at radius 2 x 10"^ G M/c? . This value has little impact 
on the result, because the main energy output of the first 
peak happens at the moment of spinar formation, and the 
spinar radius does not depend on the initial radius (see equa- 
tion 4). Simultaneously, while sampling the Scenario Ma- 
chine results, we set an upper limit on the effective Kerr 
parameter, ao ~ 44, to ensure that a core is not below the 
centrifugal barrier at the beginning of a calculation. The 
lower limit on ao is 1. Coefficient kt in the expression for 
the braking torque (5) is set to 1/3. 

Each trio of parameters, ao, Ofm, and m, yields a power 
curve with two pulses of different strength and duration. 




2 5 10 20 50 

Value of initial parameter 

Figure 2. The distribution of the pre-collapse effective Kerr pa- 
rameters (solid line) and the distribution of the masses (dashed 
Une) for the collapsing WR cores with ao > 1 for population '1' 
(see Table 1). A^total = 8000. 
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Figure 3. Energy release rate of a collapsing spinar (the sum of 
the solid and the long-dashed line). Initial parameters are m = 7, 
ao = 5, and ctm = 10~*. Time and luminosity are in the frame 
of an infinite observer. Left short-dashed profile represents char- 
acteristic observed duration of the first peak. Right short-dashed 
horizontal line is drawn at the 0.5 level from the peak value of 
Loo ■ The overlaying graph zooms in a time interval near the first 
peak. 



Various patterns can be recognized, like plateaus, tails, one- 
peak patterns, etc. 

Consider as an example the burst evolution in Fig. 3 
that shows the rate of energy release as detected by an in- 
finite observer, corrected for the gravitational redshift and 
time dilation. The solid line shows the dissipation rate of 
the kinetic energy at 'impact' (i.e., the halt at the centrifu- 
gal barrier); the long-dashed line shows the power released 
when the angular momentum is carried away. One can see 
that the two stages are clearly separated by the relative con- 
tribution of the two means of energy supply. 

The burst pattern provides two useful temporal charac- 
teristics: the separation between peaks t2 — ti and the dura- 
tion of the second peak St2. Generally, the computed dura- 
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Figure 4. Magnetic parameter versus initial effective Kerr pa- 
rameter for different classes of GRBs: with precursor (big red 
dots), with the stronger first peak (green pluses), with single peak 
(black diamonds), and merged peaks (small blue dots). Color is 
seen in the electronic version. For a comparison, the reader is 
referred to the diagram in figure 2 of Lipunov & Gorbovskoy 
(2008). 



tion of the first peak is very short, because it corresponds to 
the time-scale of a dynamic model. Thus, it should not be 
compared directly with a GRB jet duration. The process of 
a jet development or propagation is not considered or spec- 
ified here. Duration of the second peak, estimated as shown 
in Fig. 3, might also underestimate the jet duration for the 
same reason. 

We assume that a model represents a GRB with precur- 
sor if it satisfies the following condition: {Ei/5t\)/L2 < 1. 
Here Ei is the energy released during the first pulse, 
L2 = - the peak value for the second pulse. Ei is ob- 

tained by integrating over the time interval during which 
Loo is greater than I x U^^, where Z is a fractional value. 
We usually set I to 0.5 or 0.1. 

A constraint on the actual duration of the first peak 
6ti comes from observational data. It is widely accepted 
that the short and the long GRBs have different underlying 
mechanisms or progenitors (Norris et al. 2001; Balazs et al. 
2003; Fox et al. 2005). Swift data indicate that long and 
short GRBs have different redshift distributions with differ- 
ent median values of 2: 0.4 for short GRBs and 2.4 for long 
GRBs (O'Shaughnessy et al. 2008; Bagoly et al. 2006). We 
assume that both the first and the second peaks of the cal- 
culated models contribute to the class of long GRBs. Burlon 
et al. (2008) study data for Swift GRBs with precursor with 
measured redshifts, all of them having rest-frame precursor 
durations between 3 and 20 seconds. Thus, we attribute fac- 
titious durations to the first peaks, generating real values 
distributed randomly and uniformly in the above limits. 



4 RESULTS 

For each population from Table 1, we calculate 8000 models 
within wide limits of am (10~^..10~^). Among the calcu- 
lated models we can distinguish four classes of GRBs. These 
are two classes with separate peaks: GRBs with precursor 
and GRBs with a stronger first peak; third, GRBs with one 
pulse, produced by models with ao either very large or close 
to unity; the last, GRBs with very close pulses, meaning that 
the time separation between the pulses is less than the max- 
imum of two values: 0.1 s (an arbitrary chosen small value) 
and (Sti + 5t2)/2. Fig. 4 shows distributions of the initial 
model parameters for different classes of GRBs. Here precur- 
sors satisfy the following condition: 1 > {Ei/Sti)/L2 > 0.01. 
If a precursor is very weak, then a GRB is re-classified and 
becomes a single-peak GRB (diamonds in the top right cor- 
ner of Fig. 4) . 

In Fig. 5 we present characteristic time and energy dis- 
tributions for a class of models identified as GRBs with pre- 
cursor. The initial parameters for the models in Fig. 5 are 
taken from population '1' (see Sect. 2). The energies and 
durations of the peaks are calculated for a level of 0.5 of 
the peak luminosity (for illustration see Fig. 3). In Fig. 5 
the energy of pulses is given as r^iso -Bp = Eiao, where Ep 
is the calculated energy of a peak, P = 1 or 2, and Eiso is 
the rest-frame isotropic bolometric energy, which can be de- 
rived from values observed in a spectral energy band. Here 
we set r^iao = 0.01 for both peaks in order to cover with the 
modeled points the area that is occupied by the observed 
values adopted from Burlon et al. (2008) (stars in Fig. 5). 
Coefficient ?7iso incorporates our freedom in choosing level I 
as well as the actual efficiency of the gravitational collapse, 
and the fact that there are two jets. 

Burlon et al. (2008) study a sample of GRBs with pre- 
cursor activity detected by Swift and with known redshifts. 
All GRBs with precursors belong to a class of long GRBs, as 
also found for the BATSE catalogue by Koshut et al. (1995). 
We use the data from Table 1 of Burlon et al. (2008) except 
for GRB 070306 that has two precursors (see Sect. 5.1). The 
observed time intervals are recalculated to the rest-frame us- 
ing known redshifts. One can see that the separation times 
t2 — ti between the peaks calculated in the spinar model 
reproduce very well the observed time lags and that longer 
lags are also produced. From the practical point of view it is 
worth noting that separation times calculated numerically 
agree with values fam defined by (7) quite closely, within an 
accuracy of 10% for most of the models. Formula (7) can be 
rewritten as 



tan 



20- 



M 



1/3 



(8) 



10 Mq V 3 / V 10- 
In Fig. 5b, we also show by pluses results for the models 
with a stronger first peak. Some of these models can manifest 
themselves as events with a precursor due to some dispersion 
of Tjiso and differences in the efficiencies of the two peaks. The 
scatter of the observed values (stars in Fig. 5) can be also 
caused by different efficiencies in the sources. In the context 
of Fig. 5b, we notice that a correlation with slope about 1 
between emission duration time and previous quiescent time 
in the observer frame was found by Ramirez-Ruiz & Merloni 
(2001) (see their figure 3b) who investigated long and bright 
BATSE bursts containing at least one quiescent interval in 
their time history (they did not sample precursor events). 
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Figure 5. (a) Energy of the first peak vs. energy of the second peak; (b) duration of the second peak vs. separation between the 
peaks; (c) energy of precursor vs. separation; (d) precursor-to-main-peak energy ratio vs. separation. Population '1' from Table 1 is used. 
Dimensionless coefficient r)iso = 0.01 for both peaks. Magnetic parameter 10~® < Om < 10~^. Pluses (green in the electronic version) in 
(b) represent results for the models with a stronger first peak. Stars are plotted for the data taken from Burlon et al. (2008): rest- frame 
magnitudes Eiaoi mostly in 15 — 150 keV, separation times, calculated as (Ti main — Ti prcc)/(l + z), and durations of the main peaks 
.main .main )/(i + ^)- 



Their slope roughly coincides with the overal slope of the 
star-distribution in Fig. 5b. 

In Table 2 we give the approximate ratios of the num- 
ber of GRBs with precursor to the number of all GRBs for 
different binary evolution parameters. The fraction of the 
models with precursor of separation < 100 s is 8 — 10 per 
cent for most of the populations. All GRBs include two-peak 
GRBs, single peak GRBs, and GRBs with very close peaks. 
To calculate these ratios, only precursors that are 'strong 
enough to be detected' are counted; that is, they satisfy an 
arbitrary chosen condition: {E\ / St\) / L2 > 0.01. 

The numbers in the Table 2 should be regarded as rough 



estimates, as they are based on an assumption that a mod- 
eled power curve (as, for example, in Fig. 3) provides a re- 
alistic relation between Li and L2. If we halve the efficiency 
r/iso of the precursor, the relative number of precursors be- 
comes greater but not significantly (~ 10 per cent for pop- 
ulation '1', for example). To provide some glimpse of the 
significance of the numbers, we can say that, if we change 
the level for calculating peak duration and energy to / = 0.1 
and keep other parameters unchanged, then the number of 
GRBs with precursor drop to 6 per cent for population '1'. 

Normalized histograms for distributions of modeled sep- 
aration times and main (i.e., second) peak durations are 
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Table 2. Relative frequency of GRBs with precursor for different 
ranges of the modeled separation time between the precursor and 
the main peak. 



Population 



i2 - *!, s 
< 100 100 - 10^ 



100 - 10'" 



1 


0.08 


0.04 


0.10 


la 


0.11 


0.07 


0.12 


lb 


0.08 


0.09 


0.14 


2 


0.09 


0.05 


0.10 


3 


0.08 


0.05 


0.12 


4 


0.09 


0.05 


0.11 


5 


0.07 


0.05 


0.17 


6 


0.02 


<0.01 


0.03 



10 



I Modeled tg— 




10 



100 1000 10* lO"" 
Rest— frame separation, s 



10° 



Figure 6. Distribution of the separation times between the pre- 
cursor and the main pulse. Empty histograms show model esti- 
mates and are vertically scaled with the same area under the solid 
and the dashed line. Hatched histogram represents rest- frame sep- 
aration times (Ti niain ~Ti,prcc)/(l + z) from Burlon et al. (2008). 
The set of models is the same as in Fig. 5. 



plotted in Figs. 6 and 7. Observational data from Burlon 
et al. (2008) is also presented. The modeled duration of the 
second peak might be a minimum estimate of the actual 
duration because no radiation and magneto-hydrodynamic 
effects are incorporated in the present model of a GRB. 

The longest observed separation between the peaks of a 
precursor and a main pulse, ~ 400 s, was reported by Koshut 
et al. (1995) but no conclusion about the rest-frame time 
separation can be made without knowing z. Nevertheless, 
our results are in apparent agreement with the data from 
figure 13 of Koshut et al. (1995) that shows the ratio of the 
precursor total counts to the main episode total counts ver- 
sus separation time between detectable emission. The dis- 
tribution in Fig. 8 is constructed for population '1' for anal- 
ogous values where we try to emulate cosmologically dilated 
separation times by multiplying all values {t2 — t\) by a con- 
stant factor. 

Summarizing, the spinar paradigm explains observed 
precursor separation times and predicts that there are very 
early precursors with rest-frame time separations more than 
100 s, which is the maximum determined rest-frame separa- 
tion so far (Burlon et al. 2008). Our results show that the 
rate of GRBs with precursor, which have separations longer 



10 - 




10 100 
Rest-frame duration, s 



1000 



Figure 7. Distribution of the duration of the main pulse. Empty 
histogram shows model results for events with precursors and is 
vertically scaled. Hatched histogram represents rest-frame dura- 
tion times (T2 main — Ti main)/(l + ^) taken from Burlon et al. 
(2008). The set of models is the same as in Fig. 5. 
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200 
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Figure 8. Ratio of the energy released during the precursor and 
the main peak versus separation time multiplied by the constant 
factor. Dots represent model results for GRBs with precursor for 
population T'. The average value of the redshift for long GRBs 
ziang = 2.4. Stars are adopted from figure 13 of Koshut et al. 
(1995) that shows the ratio of the precursor total counts to the 
main episode total counts versus separation time between de- 
tectable emission for BATSE GRBs with precursor. 



by up to one order of magnitude than those observed so 
far, is about half the detected rate (compare the second and 
third columns of Table 2). The number of modeled GRBs 
with precursor with lead times in the range from 100 s to 
10^ s is comparable to the frequency of precursor events ob- 
served so far. At the same time, the expected average energy 
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output of a precursor drops with increasing separation time 
as illustrated by Fig. 5c. 

One should keep in mind that a rest-frame time inter- 
val as, for example, 100 s, is dilated to {1 + z) x 100 s in 
the observer's frame. At this point, an important general 
conclusion can be made: GRBs at high redshifts will show 
precursor events which should arrive {1 + z) x 100 ~ 1000 s 
in advance of the main event. This prediction is very strong 
as it does not depend on a GRB model. 



5 DISCUSSION 

The relative frequencies of the GRBs with precursor, ob- 
tained in the spinar paradigm, agree reasonably well with 
observational facts. Koshut et al. (1995) find precursors in 
~ 3 per cent of BATSElong and short GRBs detected before 
May 1994; Lazzati (2005) finds that precursors are present 
in ~ 20 per cent of bright, long GRBs in the final BATSE 
catalog; Burlon et al. (2008) obtains that about 14 per cent 
from 105 GRBs with measured redshift, observed by Swift 
before March 2008, have precursor activity. Note that these 
researches use different criteria for precursor search and se- 
lection. All these data agree quite well with the modeled ra- 
tios within an accuracy of 5 — 10 per cent (see Table 2). Note 
that there is no significant diff'erence in the ratios between 
the different binary evolution scenarios (for wind model Wl) 
and different ki. Apparently, the strong wind scenario expe- 
riences the most difficulty in describing the reported ratios 
(population '6'). 

To speak of a more robust result, the spinar model pro- 
vides a clear explanation of the diversity of observed time 
separtions between precursor and main peak, which, more- 
over, can be related to the main physical parameters of a 
collapsing object in a simple way (equation 7 or 8). To our 
knowledge, explaining both relatively short (~ 10 s) and 
long (~ 100 s) switch-off times presented difficulties for most 
of the GRB engine models (see Wang & Moszaros 2007; 
Drago & Pagliara 2007, and references therein). 

Accretion has so far been neglected by us, but could 
be easily incorporated in the equations, as pointed out 
by Lipunov & Gorbovskoy (2008); this is a topic for a forth- 
coming paper. As a preliminary result, we found that ac- 
cretion of up to > 1 M.Q during the overall GRB time 
(10 — 100 s) does not change qualitatively any of the results 
presented. 

5.1 Multi-precursors 

Among 15 GRBs with precursor studied by Burlon et al. 
(2008) there is one (GRB 070306) with two precursors. The 
first precursor is ~ 105 s from the second precursor, which 
is ahead of the main pulse by ~ 95 s. Each pulse is stronger 
than the previous one. This GRB is not shown in the figures 
in the present paper. 

One way to explain such a phenomenon is an unstable 
regime of jet with a constantly operating central engine. 
A mechanism of a modulated relativistic wind, resulting in 
long time gaps in emission, was considered by Ramirez-Ruiz 
et al. (2001). 

Other possibilities appear if we consider the basic el- 
ements of the spinar paradigm. As we argue in the Intro- 



duction, the course of a collapse is affected by the initial 
distribution of the angular momentum in the rotating ob- 
ject. In this regard, spinar and coUapsar models are 'oppo- 
site' options. In reality there is most probably an assortment 
of mixed scenarios. One possibility, provided by a particu- 
lar angular momentum distribution in a pre-coUapse core, is 
that a black hole forms first (primary pulse), then a spinar 
(heavy accretion disc) forms around it (second pulse), fol- 
lowed by a fatal collapse (hyper-accretion) of a spinar into 
a black hole (third pulse). Another explanation for a throe- 
stage collapse can be that an accretion disc aggregates at 
some stage of a spinar evolution. Here we should mention 
the 'initial jet and fallback collapsar' scenario (Wang & 
Mcszaros 2007) with its characteristic disc time of ~ 100 s. 
In its context, we can speculate why the majority of GRBs 
with precursor (or GRBs with one or more quiescent peri- 
ods, see Ramirez-Ruiz & Merloni (2001)) are not three-step 
events. The fallback scenario works for explosions that are 
not too weak or too strong according to Fryer (1999). An ex- 
plosion which is too weak leads to a direct collapse to a BH, 
or a spinar if it rotates fast enough. Rotating progenitors 
possibly favor 'direct collapse to a spinar' without a huge 
explosion; this view is supported by general and computa- 
tional (Monchmeyor 1991; Yamada & Sato 1994) arguments 
that rotation weakens the bounce and hence the explosion. 

Numerical studies yield further options for gravitational 
collapse: fragmentation of a rotating collapsing body, a mul- 
tiple black hole system formation (sec, for example, Berezin- 
skii et al. 1988; Imshennik 1992; Zink et al. 2007), and its 
eventual merger. Such scenarios (see also King et al. 2005) 
can provide a burst with several pulses, but are beyond the 
scope of the present work. 

5.2 Lower limit on a precursor energy 

Let us provide some relations for the energetics of the pre- 
cursor jet. The precursor jet is apparently the first to break 
through a shell that probably surrounds the central object. 
Lipunov k, Gorbovskoy (2008) arrive at a condition for a jet 
to penetrate the surrounding shell. The jet breaks through 
if the momentum imparted on a part of the shell is greater 
than the momentum corresponding to the escape velocity: 

o 7-1 7-,p ^ , , / 2GAicore /„N 

where 13 is the efficiency of processing the energy of the 
collapse to the pushing jet, c is the speed of the jet, which 
is approximately the speed of light, Os/47r is the portion of 
the shell's surface subject to the jet. Substituting universal 
constants, one has: 

0:011^™ (,27To^j •^'"S' 

(10) 

where msheii is in solar masses. 

Now consider the relations for the observed jet: 

■E^jS = ^7 -S'j^ti ^iet = -Biso ■ (11) 

It is natural to assume that the energy of the jet punching 
the shell is not the same as that of the observed 7-jet and 
to introduce a factor fc-^ < 1 between them, which depends, 
among other things, on the spectral band observed. Here 
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Qob is the solid angle corresponding to the observed jet's 
opening angle, Eiso is the rest-frame isotropic bolometric 
energy, which can be derived from the values observed in a 
spectral energy band. 

From the above wo derive the minimum value Eiso that 
corresponds to a jet still capable of going through the shell: 

Eiao{min) = 6 x Kf^k-, msheii ^ 2 >f los'r ) 

(12) 

Let us now illustrate this with some numbers. Fig. 5a 
possibly indicates that Eisoi'inin) ~ 10^" erg from the data 
of Burlon et al. (2008). Comparing with equation (12), one 
can obtain an implicit estimate of k-y for the observational 
data mentioned as low as 0.1 — 1 per cent, if Hs ~ f^ob and 

mshell ~ 1. 

In Sect. 4, we introduced the parameter r^iao to match 
the modeled energy of the pulses with the rest-frame isotrop- 
ical values Eiso- We estimated it as r?iso ~ 0.01. Using (9) and 
(11), rjiso can be expressed as follows: 

_ Eiso 47r . , , 

If the half-opening angle of a jet ~ 0.1 rad (e.g., Meszaros 
2006) and k-y ~ 0.001 - 0.01 then ilob ~ 0.03 sr ~ 47r/400 
and we come at an estimate for the effleciency of processing 
the energy of the collapse to the jet in the shell: /? ~ rjiso ~ 
0.01. 



6 CONCLUSION 

Currently, much effort is put into numerical computations 
of GRBs and simulations of jets emerging from rotating 
magnetized configurations, to name just a few, by Lyutikov 
(2006); Barkov & Komissarov (2008); McKinney & Bland- 
ford (2009); Takiwaki ct al. (2009) and others. We believe 
that the progress with the spinar mechanism for GRBs de- 
mands a further development by means of a MHD simular 
tion. 

In the present work, we calculate many models of two- 
stage spinar evolution provided that fast rotating stellar 
cores are produced in binary systems due to a tidal synchro- 
nization mechanism. The pre-collapse parameters of WR 
stars are calculated by the Scenario Machine dedicated to 
the population synthesis of binary systems (Lipunov et al. 
1996a, 2007). 

We analyse the resulting set of GRBs, which can be 
classified into different classes. For GRBs with precursor, 
separation times between the precursor and the main pulse 
are fully consistent with observational data. We find that 
precursor events with the rest-frame time separation from 
the main pulse < 100 s occur in about 10 per cent of the 
modeled GRB events, for the weak wind scenarios of bi- 
nary evolution and independently of the poorly established 
value of kj. This rate, however, depends on the specifics of 
the jet generation and its reliability is limited by the frame- 
work of our simple model. Taking this into account, we think 
the modeled rate agrees very well with the observed values. 
Comparison of the time and energetic characteristics of the 
modeled bursts with the observational data for GRBs with 
measured redshifts (Burlon et al. 2008) is made and a good 
general agreement is found. 



We predict that very early precursors with a rest-frame 
time separation more than 100 s should exist, and the num- 
ber of GRBs with precursor with separations in the range 
from 100 s to 1000 s is about half the number of those ob- 
served so far. Super-early primary pulses, up to 10® s, are 
also found in the model; they, however, may be too weak to 
be detected or to develop as a jet. 

Finally, whether or not the present model works in all 
cases, the observational evidence is that rest-frame separa- 
tion times of ~ 100 s take place. Higli-rcdsliift long GRBs, 
up to « > 10, are expected to be detected by future exper- 
iments (e.g., Salvaterra et al. 2008). Thus, we confidently 
expect precursor events arriving up to >1000 seconds in ad- 
vance of the main GRB episode. 
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